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C^ I We present a new continuous data assimilation algorithm based on 

ideas that have been developed for designing finite-dimensional feed- 
back controls for dissipative dynamical systems, in particular, in the 
context of the incompressible two-dimensional Navier-Stokes equa- 
^ ■ tions. These ideas are motivated by the fact that dissipative dynami- 

Qv . cal systems possess finite numbers of determining parameters (degrees 

Q^ I of freedom) such as modes, nodes and local spatial averages which gov- 

^— ^ ' ern their long-term behavior. Therefore, our algorithm allows the use 

'nJ" . of any type of measurement data for which a general type of approxi- 

Zfi '. mation interpolation operator exists. Our main result provides condi- 

tions, on the finite-dimensional spatial resolution of the collected data, 
sufficient to guarantee that the approximating solution, obtained by 
our algorithm from the measurement data, converges to the unknown 
reference solution over time. Our algorithm is also applicable in the 
C^_' context of signal synchronization in which one can recover, asymptoti- 

cally in time, the solution (signal) of the underlying dissipative system 
that is corresponding to a continuously transmitted partial data. 
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1 Introduction 

The goal of continuous data assimilation, and signal synchronization, is to 
use low spatial resolution observational measurements, obtained continuously 
in time, to accurately find the corresponding reference solution from which 
future predictions can be made. The motivating application of continuous 
data assimilation is weather prediction. The classical method of continuous 
data assimilation, see Daley [TU], is to insert observational measurements di- 
rectly into a model as the latter is being integrated in time. We propose a 
new approach based on ideas from control theory, see Azouani and Titi [2]. 
A slightly similar approach in the context of stochastic differential equations, 
using the low Fourier modes as observables/measurements, appears in a re- 
cent work by Blomker, Law, Stuart and Zygalakis [1] . Rather than inserting 
the measurements directly into the model, i.e. into the nonlinear term, we 
introduce a feedback control term that forces the model toward the reference 
solution that is corresponding to the observations. This is motivated by the 
fact that the measured data is usually obtained as the values of the exact 
solutions at a discrete set of spatial nodal points, and that it difficult to insert 
this data directly into the underlying equation, because it is not possible to 
obtain the exact values of the spatial derivatives. One should observe that 
in order to guarantee a unique corresponding reference solution one has to 
supply observational data with enough spatial resolution. This is the object 
of this paper. 

While the classical method of continuous data assimilation is simple in 
concept, special care has to be taken concerning how the observations are 
inserted into a model in practice. For example, it is generally necessary to 
separate the fast and slow parts of a solution before inserting the observations 
into the model. The method proposed here does not require such a decom- 
position. Since the observations are not directly inserted into the model, we 
can rely on the dissipation already present in the dynamics to filter the ob- 
served data, i.e. the viscous term will suppress the "spill over" oscillations in 
the fine scales. The advantage of this approach is that it works for a general 



class of interpolant observables without modification. 

Let u{t) represent tlie state at time t of the dynamical system in which 
we are interested, and let Ih{u{t)) represent our observations of this system 
at a coarse spatial resolution of size h. Given observational measurements, 
Ih{u{t)), for t G [0, T], our goal is to construct an increasingly accurate initial 
condition from which predictions of u{t), for t > T, can be made. We do 
this by constructing an approximate solution v{t) that converges to u{t) over 
time. 

Suppose the time evolution of u is governed by a given evolution equation 
of the form 

where the initial data, Uq, is missing. Our algorithm for constructing v{t) 
from the observational measurements Ih{u{t)) for t G [0,T] is given by 

'^^=F{v)~^^h{v)+^^h{u), (2) 

^(0) = i;o, (3) 

where /x is a positive relaxation parameter, which relaxes the coarse spatial 
scales of v toward those of the observed data, and vq is taken to be arbi- 
trary. It is worth stressing that our algorithm is designed to work for general 
dissipative dynamical systems of the form ([T]). Such systems are known to 
have finitely many degrees of freedom in the form of determining parame- 
ters of the type Ih{u), see, for example, Cockburn, Jones and Titi [8], Foias, 
Manley, Rosa and Temam [13], Foias and Prodi [H], Foias and Temam [T5] . 
fT6], Jones and Titi [21], [22], and references therein. The incompressible 
two-dimensional Navier-Stokes equations provide a concrete example of a 
dissipative dynamical system of this type. 

We consider here the incompressible two-dimensional Navier-Stokes equa- 
tions, as a paradigm, because they are amenable to mathematical analysis 
while at the same time similar to the equations used in realistic weather 
models. Thus, we shall suppose the evolution of u is is governed by the 
Navier-Stokes system 

Oil 

— - uAu + (n • V)n + Vp = f (4) 

V ■ n = 0, (5) 



in the physical domain Q, with either no-shp Dirichlet, or periodic, boundary 
conditions. Here u{x,t) represents velocity of the fluid at time t at position 
X, p > represents the kinematic viscosity, p{x,t) is the pressure and f{x,t) 
is a time dependent body force applied to the fluid. 

In the case of no-slip Dirichlet boundary conditions we take u = on 
(9^2. The domain fi is an open, bounded and connected set in R^ with C^ 
boundary, such that dQ can be represented locally as the graph of a C^ 
function. In the case of periodic boundary conditions we require u and / to 
be L— periodic, in both x and y directions, and take Q = [0,L]^ to be the 
fundamental periodic domain. 

Continuous data assimilation, in the context of the incompressible two- 
dimensional Navier-Stokes equations, was flrst studied by Browning, Hen- 
shaw and Kreiss in [7], later by Henshaw, Kreiss and Ystrom in [20] and 
also by Olson and Titi in |23] and [25j, motivated by the concept of flnite 
number of determining modes which was introduced for the flrst time in 
[H] (see also [I3], |2l], and references therein). These studies treated the 
case of periodic boundary conditions, where the observations were given by 
the low Fourier modes with wave numbers k, such that \k\ < 1/h. Since 
the low modes essentially represent the large spatial scales of the solution, 
the classical data assimilation algorithm works well for this type of obser- 
vations. In addition, it is worth mentioning that this method is consistent 
with some of the signal synchronization algorithms. Most recently, similar 
idea has also been introduced in [12] to show that the long-time dynam- 
ics of the two-dimensional Navier-Stokes equations can be imbedded in an 
inflnite-dimensional dynamical system that is induced by an ordinary differ- 
ential equations, named determining form, which is governed by a globally 
Lipschitz vector fleld. 

The method of constructing v, given by ([2]), allows the use of general 
interpolant observables, given by interpolants Ih'-H^{Q) — )■ L'^{Q) that are 
linear and satisfy the following approximation property: 

11^ - h{<^)\\l^n) < coh^lMmin) (6) 

for every (p G H^{Q). The orthogonal projection onto the low Fourier modes, 
with wave numbers k such that \k\ < 1/h, mentioned above, is an example 
of such interpolant observable. However, there are many other interpolant 
observables which satisfy ([6]). 

One physically relevant example of an interpolant which satisfles condi- 
tion (El) are the volume elements studied in [211 and 1221 (see also Foias and 



Titi [T7]). In this case 

Ih{v{x)) = Y.ifjXQ,{x) where Vj = j^ J ^{x)dx, 



^ N 



j=i 



and the domain fl = [0,L]^, for the periodic boundary conditions case, has 
been divided into A^ equal squares Qj, with sides h = L/\/N. Volume 
elements generalize to any domain Q on which the Bramble-Hilbert lemma 
holds. An elementary discussion of this lemma in the context of finite element 
methods appears in Brenner and Scott [5J. 

In addition, we also consider interpolant observables given by linear inter- 
polants Ih'-H'^i^) -> L^{Q), that satisfy the following approximation prop- 
erty: 

ll<^ - 4(V^)||i2(n) < ^C^/l^||^||^2(f^), (7) 

for every cp G H'^{Q). An example of this type of interpolant is given by 
measurements at a discrete set of nodal points in Q. Specifically, let h > 
be given, and let Q, = Uj^iQj, where Qj are disjoint subsets such that 
diamQj < h, for j = 1, 2, . . . , Nh, and let Xj G Qj be arbitrary points. Then 
set, for example, 

hi^ix)) = J2 fixk)XQjix). (8) 

fc=l 

Following ideas in [22] (see also [16]) one can show that Ihiy^) satisfies ([7]). 

Our paper is organized as follows. First, we recall the functional setting 
of the two-dimensional Navier-Stokes equations necessary for our analysis 
and then use this setting to formulate our new method of continuous data 
assimilation. After this we proceed to our main task, that of finding con- 
ditions under which the approximate solution obtained by this algorithm of 
data assimilation converges to the reference solution over time. Section 3 
treats the case of smooth, bounded domains with no-slip Dirichlet boundary 
conditions, while section 4 treats the case of periodic boundary conditions. 
Our main results may be stated as follows: 

Theorem 1. Let Q be an open, bounded and connected set in R^ with C^ 
boundary, and let u be a solution to equations 0)-([^ with no-slip Dirichlet 
boundary conditions. Assume that I^ satisfies ^, with h small enough such 
that 

1/h^ > ciXiG^, 



where C\ is a constant given in [3&\) . Then there exists fi > 0, given explicitly 
in PropositionUl such that \\v — u\\]^2(^q^ — )■ exponentially, as t —^ oo. 

Here G denotes the Grashof number 

G=^limsup||/(t)|U2 (9) 

where Ai is the smallest eigenvalue of the Stokes operator subject to homoge- 
neous Dirichlet boundary conditions. Let us remark, again, that the constant 
Ci is depends only on cq, given in IQ, and the shape, but not the size, of the 
domain Q. In particular, ci is given by ([3SD where the constant c is chosen 
so the bound f lTB]l on the non-linear term holds. Moreover, /i may be chosen 
equal to bc^G^uXi as indicated in Proposition [H below. 

Results similar to Theorem [2] hold when I^ satisfies ([7]), however, we omit 
the proof of this result in the case of no-slip Dirichlet boundary conditions 
and instead proceed directly to the the case of periodic boundary conditions 
where sharper estimates may be obtained. In particular, we prove 

Theorem 2. Let Q = [0, L]^ and let u be a solution to equations ^-^ with 
periodic boundary conditions. Let I^ satisfy either ^ or ^, with h small 
enough such that 

l/h^>C2XiGil + \ogil + G)), 

where C2 is a constant given in / T?^) . Then there exists fi > 0, given explicitly 
in Proposition^^ such that \\v — u\\h^{q) — ?■ exponentially, as t -^ oo. 

Let us remark again that C2 depends only on cq, and that /i may be chosen 
as 3c2i'\iG{l + log(l + G))/co. In particular, /i is given in Proposition |2] and 
02 is defined in (1361) as an increasing function of cq and c, where c is chosen 
large enough so that the bounds in both fl2^ and flMl) hold. 

Note that the estimate on the length scale h in Theorem [2] can be com- 
pared to previous results reported in [23]. Let v(t) be the approximate solu- 
tion obtained by the method of continuous data assimilation introduced in 
[23] for the interpolant observable Ih{u) given by projection onto the Fourier 
modes with wave numbers \k\ < 1/h. In [24J it was shown, that for small 
values of h, such that l//i^ ~ AiG, \\u{t) — 'y(t)||ffi{o) -^ exponentially fast, 
as t — )■ oo. Up to a logarithmic correction term. Theorem |2] states similar 
estimates on h for the new algorithm which covers a much wider class of 
interpolant observables. 

The final section of this paper discusses numerical simulations, which are 
in progress, related works, and closes with a few concluding remarks. 



2 Preliminaries 

This section reviews the functional setting of the two-dimensional Navier- 
Stokes equations with no-slip and periodic boundary conditions, recalls some 
facts that will be used in the remainder of the paper and then gives an 
explicit formulation of our new method for continuous data assimilation in 
this context. Following Constantin and Foias |9], Foias, Manley, Rosa and 
Temam [13j, Robinson [26j and Temam [27j, we begin by defining a suit- 
able domain Q and space V of smooth functions which satisfy each type of 
boundary conditions. 

No-slip Dirichlet Boundary Conditions. Let Q be an open, bounded 
and connected domain with C^ boundary. Define V to be set of all C°° 
compactly supported vector fields from Q to R^ that are divergence free. 

Periodic Boundary Conditions. Let Q = [0, L]^ for some fixed L > 0. 
Define V to be the set of all L-periodic trigonometric polynomials from R^ 
to R^ that are divergence free and have zero averages. 

Given V corresponding to either type of boundary conditions let H be 
the closure of V in L^{n; R^) and V be the closure of V in H^{n; R^). The 
spaces H and V are Hilbert spaces with inner products 

(u,v) = / u(x)-v(x)dx and ((u,v)) = T^ / —^—-^dx, 
Jn f-^^ Jn oxj oxj 

respectively. Denote the norms of H and V by 

\u\ = \J{u,u) and ||u|| = ^{{u,u)), 

and the dual of V by V* with the pairing {u, v) where u &V* and v &V . 

Define the Leray projector Po- as the orthogonal projection from L^(r2; R^) 
onto H ^ and define the Stokes operator A:V ^ V* , and the bilinear term 
B:V ^V ^V* iohe the continuous extensions of the operators given by 

Au = -Pa^u and B{u, v) = P„{u ■ Vv), 

respectively, for any smooth solenoidal vector fields u and v in V. 

Denote the domain of A by T>{A) = {u ^ V : Au E H}. The linear oper- 
ator A is self-adjoint and positive definite with compact inverse A~^: H ^ H. 
Thus, there exists a complete orthonormal set of eigenfunction Wi in H such 



that Awi = XiWi where < Aj < Aj+i for i G N. Writing Ai as the smallest 
eigenvalue of A we have the following Poincare inequalities: 

a u e V then Xi\u\^ < \\uf , (10) 

if M e V{A) then \i\\uf < |AMp. (11) 

Note that for u G H, \u\ = ||^^||L2(f7) and for u G V the Poincare inequality 
implies ||m|| is equivalent to ||u||_f/i(Q). 

The bilinear term B has the algebraic property that 

{B{u,v),w) = -{B{u,w),v) (12) 

for u,v,w E V, and consequently the orthogonality property that 

{B{u,w),w) = 0. (13) 

Here the pairing (•, •) denotes the dual action of V* on V. Details may be 
found, e.g., in [9], [I3], [26] and [27]. 

In the case of periodic boundary conditions the bilinear term possesses 
the additional orthogonality property 

{B{w,w),Aw) = 0, for every weV{A); (14) 

and consequently one has 

{B{u, w), Aw) + {B{w, u), Aw) = -{B{w, w), Au). (15) 

Note that the bilinear term satisfies a number of inequalities which hold 
for either no-slip or periodic boundary conditions. These are 

\{B{u,v),w)\ < c\u\^/^\\u\\^^^\\v\\\w\'^''\\w\\^/\ (16) 

for every u,v,w G V, 

\{B{U,V),W)\ < c|m|1/2||^||1/2||^||1/2|^^|1/2|^| (^7^ 

for every u E V, v E 1^{A) and w E H, and 

\{B{u,v),w)\ < c\u\^^^\Au\^/^\\v\\\w\, (18) 

for every u G ^(A) and v,w G V, where c is a dimensionless constant 
depending only on the shape, but not the size, of Q. These inequalities 
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may be obtained from the Holder's inequality, the Sobolev inequalities and 
Ladyzhenskaya's inequality, see, e.g., |H], [13], [2S] and [27]. 

With the above notation we write the incompressible two-dimensional 
Navier-Stokes equations in functional form as 

— + uAu + B{u,u) = f (19) 

with initial condition u{0) = uq. We have assumed f E H so that P^f = f. 
As shown in [9], [I3], [26] and [27] these equations are well-posed; and possess 
a compact finite-dimensional global attractor, when / is time-independent. 
Specifically, we have 

Theorem 3 (Existence and Uniqueness of Strong Solutions). Suppose uq G 
V and f G L°°{{0, oo), H). Then the initial value problem /[T^) has a unique 
solution that satisfies 

fill 

ueC{[0,T];V)nL\iO,Ty,DiA)) and - e L\{<d,T)-H\ 

for any T > 0. 

We now give bounds on solutions u of 0191) that will be used in our later 
analysis. With the exception of inequality fl22]) due to Dascaliuc, Foias and 
Jolly [11] these estimates appear in any the references listed above. 

Theorem 4. Fix T > 0, and let G be the Grashof number given in (0). 
Suppose that u is the solution of / fiPj) . corresponding to the initial value uq, 
then there exists a time to, which depends on uo, such that for all t > to we 
have: 



rt+T 

\u{t)\'^ <2u^G^ and / ||M(r)f rfr < 2(1 + Tz/AOz/G^. (20) 



h(t)f <2^/^AlG'^ [^ \Au{T)\^dT<2{l + TuXi)uXiG^; (21) 



rt+T 
It 
In the case of periodic boundary conditions we also have: 

rt+T 
It 
furthermore, if f & H is time-independent then 

\Au{t)\^<cu^Xl{l + G)\ (22) 



We now write the continuous data assimilation equations ([2]) for the in- 
compressible two-dimensional Navier-Stokes equations. Let m be a strong 
solution of ©-Q, or equivalently flT9|) . as given by Theorem [3l and let Ih 
be an interpolation operator satisfying ([6]) or ([7]). Suppose that u is to be 
recovered from the observational measurements Ih{u{t)), that have been con- 
tinuously recorded for times t in [0,T]. Then, the approximating solution v 
with initial condition vq G V, chosen arbitrarily, shall be given by 

dv 

— -iyAv + {vV)v + Vq = f + fiihiu) - h{v)), 

on the interval [0, T]. Using the above functional setting the above system is 
equivalent to 

— + iyAv + B{v, v) = f + fiPMhiu) - h{v)), (23) 

on the interval [0,T]. 

If we knew uq exactly, then we could take vq = uq and the resulting 
solution V would be identical to u for all time; this is due to the uniqueness of 
the solutions of fl23l) (see Theorem [5|, below). However, if we knew Uq exactly, 
there would be no need for continuous data assimilation in the first place and 
one could integrate (IT^ directly with the initial value Uq. Intuitively speaking 
it makes sense to take Vq = Pfjlh{u{0)), which is the initial observation of 
the solution u. However, vq chosen in this way may not be an element 
of V. The main point of the data assimilation method given in (123|) is to 
avoid the difficulties which come from the direct insertion of observational 
measurements into the approximate solution. A choice for vq in agreement 
with this philosophy is Wq = 0. In fact, our results to hold equally well when vq 
is chosen to be any element of V. In either case we obtain an approximating 
solution V constructed using only the observations of the solution Ih{u) and 
the known values of z/ and /. 

We now show the data assimilation equations fl25]) are well-posed. When 
Ih satisfies ([H]) we show well-posedness for both no-slip Dirichlet and periodic 
boundary conditions. When Ih satisfies ([7]) we will deal here, for simplicity, 
with only the case of periodic boundary conditions. 
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Theorem 5. Suppose Ih satisfies ^ and fiCoh"^ < v, where Cq is the con- 
stant appearing in ^. Then the continuous data assimilation equations /[^) 
possess unique strong solutions that satisfy 

nil 

veCi[0,T];V)nL'iiO,Ty,DiA)) and - e L\{Q,T)-H), (24) 

for any T > 0. Furthermore, this solution depends continuously on the initial 
data vq in the V norm. 

Proof. Define g = f + ^Palhiu). Theorem [3] implies u G C([0,T]; V). Con- 
sequently 

\Pah{u)\ <\U- h{u)\ + \U\ < {c'J'h + X^'^')\\U\\ 

implies that P^4(m) G C{[0,T];H). Hence g G C{[0,T];H). This means 
there is a constant M such that {gl"^ < M for every t G [0,T]. 

We now show the existence of solutions v to (!23|) using the Galerkin 
method. The proof follows the same ideas as the proof of Theorem [31 Let Pn 
be the n-th Galerkin projector and f " be the solution to the finite-dimensional 
Galerkin truncation 

dv"- 

+ j,Av"+PnB{v'',Vn=Png-fiPnIh{vn ^ ^ 

dt (25) 

i;"(0) = PnVo. 



First, we observe that ( !25l) is a finite system of ODEs, which has short time 
existence and uniqueness. We focus on the maximal interval of existence, 
[0,T„), and show uniform bound for Vn, which are independent of n. This 
in turn will imply the global existence for (!25|) . Thus, our aim is to find 
bounds on f " which are uniform in n. This will then show global existence 
of solutions to (!23|) . In the estimates that follow, we denote the Galerkin 
solution f " by v for notational simplicity. 

Begin by taking inner products of (1251) with v to obtain 

2^I^P + '^ll^f = {9,v) - fi{Ih{v),v) 

= {g, v) + /i(f - Ihiv),v) - /i|t;p 

< — \q\ H \\v\\ . 

- 2/i'^' 2 " " 

11 



By hypothesis h is so small that ficoh^ < v. Therefore, 



— \v\'^ + v\\vf<-\g\^, (26) 

at jji 



and consequently 



— |wp + v\i\v\^ < -M, for every t G [0,T„). (27) 

Multiplying ( 1271) by e^^^^ and integrating yields 



2 , ,2 ^ 



where 

As this bound holds uniformly in n for T„ arbitrarily large, we have global 
existence on the interval [0,T], for all T > 0. Now, integrating (12B]) yields 

t 
It follows that 



I^WI'-l^ol' + ^^Tll^f <-M. 
JO u 



/ ||f(r)||(ir < ay, for every t G [0,T], 
Jo 



where 

1 T 

^y = -|t'oP + — M. 

z/ /iZ/ 

Now, take inner products of (^^ with Av to obtain 

--||T;f + ulAv]"" + {B{v, v), Av) = {g, Av) - /i(4(t;), Av). 



Inequality ( !T7j) implies 

|(5(t;,t;),Ai;)| < c|t;|^/lt;|||Ai;|=^/2 



1 /63/4 , .,„, uV 3 /z/3/4 y/3 



< ^ K^ K^ + o Mit^ • 
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Furthermore, 

\{g,Av)\<\g\\Av\<-\g\^+'^\Av\^ 



V 



and by ([6]) along with the assumption that [ic^b? < i^ we obtain 



7/ /I 



V ' ' " 4' 

/A^ 2 , '^M |2 II ||2 ^^\ A |2 

< t> H — At; — u k; < —\Av\ . 

- ^ II II ^11 /-ii II - 41 I 



Therefore, 

(it" " ' ' ' — 2/3 I I II II ' ^ 



^ii^^f +H^^i'<^i^ni^r+-i^i', (28) 



and consequently 

^ h^ ^ K^ K^ < -^ < -^, 29 

at v^ V V 

for every t G [0,T]. Define 

V^(t)=exp|-^|jt;n|t;f|. (30) 

Since 

I \ \2\\ \\2 ^ 2 I w \\2 ^ 2 2 ^ £ 4- ^ \r\ '-r^ 

/ 1^1 ll^ll ^ Ph II"^!! ^ PtfCTy < 00, tor every rG[U,iJ, 

JO JO 

we have that ipit) > for every t G [0,T]. Multiplying ( |27|) by '?/'(t) and 
integrating yields 



\\v{t)r< 

where 



<-7^{lkoir + -M^ ^(s)rfs|<p^, for all tG[0,T], 



2 1 fn 1,2 4T,^ 



Now, integrating (!28|) yields 



JQ I/'^ Jo ^ V ^ 



.2 
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for every t G [0, T], where 

54c^T 



^nA) - —jj- 



p^pl^ + ^My 



The bounds pv and o"x)(a) are uniform in n. Uniform bounds on \dv/dt\ 
then proceed in exactly the same way as for the two-dimensional Navier- 
Stokes equations. Since the estimates on the Galerkin solutions are uniform 
in n, Aubin's compactness theorem pQ allows one to extract subsequences in 
such a way that the limit v satisfies (123|) and i^^. 

Next, we show that such solutions are unique and depend continuously 
on the initial data. Let vi and V2 be two solutions for ( 123|) both satisfying 
the conditions in f l24p . Choose K large enough such that ||fi|P < K and 
||'W2|P < K for almost every t G [0,T]. Let 6 = Vi — V2- Then 5 satisfies 

— + PA5 + B{v^,6) + B{6, V2) = -iJiPMS). 
at 

Taking inner product with A5 yields 

ilpf + ^\A5\^ + ^B{v,, 6), A6) + {B{6, V2), A6) = -p{h{S), A6). 
Here we used the fact that 

which can be justified by Lemma 1.2 in Chapter 3 of Temam [28] or Theorem 
7.2 in Robinson [SB] which is due to Lions-Magenes [23] ■ Estimate the right- 
hand side of this equation as 

~pihi5), AS) = p{5 - h{5),A5) - p\\5f 

<^l^\\Sf + ^_\A6\^^p\\6r<''-\A6\\ 
where we have again used the hypothesis that pc^h? < v. It follows that 

^^ll^f + ^l^5p<l(i?(t'i,5),^5)l + l(5(5,^2),A5)|. (31) 
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The proof of uniqueness and continuity now proceeds as for the two- 
dimensional Navier-Stokes equations. In particular, estimate the non-linear 
terms on the right-hand side of fl3T]) using flTTl) and flTSj) as 



and 



\{B{5^v,lA5)\ < c\5\'''\\v,\\\A5f' < ^J^Pf + -\A5\\ 



Therefore, 



l"^"'-^^"^"'' ^°'^" '^[°'^]- 



Integrating yields 



2 „,„2 r27c^K2 



Thus, the solutions v to ( !23l) . which satisfy (!24l) . also satisfy i; G C([0, T], \/), 
and depend continuously on the initial data in the V norm. D 

Theorem 6. In the case of periodic boundary conditions suppose that Ih sat- 
isfies ^, and ficoh^ < v, where Cq is the constant appearing in ^. Then the 
continuous data assimilation equations / f^) possess unique strong solutions 
that satisfy [24\ ), for any T > 0. Furthermore, this solution is in C{[0,T], V) 
and depends continuously on the initial data vq in the V norm. 

Proof. The proof is similar to the proof of Theorem \5\ but makes use of the 
identity f fT4|) to obtain estimates on ||i;|| and /q \Av\'^ directly. D 

The algorithm given by equation (!23l) for constructing the approximate 
solution V contains two parameters h and /i. The first parameter h has di- 
mensions of length and corresponds to the resolution of the observational 
measurements represented by Ih{u). Smaller values of h correspond to spa- 
tially more accurate resolved measurements. The relaxation parameter ix 
controls the rate at which the approximating solution v is forced toward the 
observable part of the reference solution u. Larger values of /x cause Ih{y) 
to faster track Ih{u). It is the parameter /i which distinguishes fl23|) from 
the previous methods of continuous data assimilation studied in [7], [20] 
and [2g. 
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The condition that ficoh"^ < u, given in Theorem El places a restriction 
on the size of fih"^ compared to the viscosity u, sufficient to ensure the data 
assimilation equations are well-posed. This restriction is due to the fact 
that the the interpolant operator filh might generate large gradients and 
spatial oscillations ( "spill over" to the fine scales) that need to be controlled 
(suppressed) by the viscosity term. Notice that in the case when Ih = Pnih-i 
where P^h is the orthogonal projection onto the linear sub-space spanned 
by the Fourier modes with wave numbers \k\ < rrih = j^, such oscillations 
are not generated, since —{fj,Ih{v),v) = —^iPm^vl"^ and —{fiIh{v),Av) = 
—liWPm^vW^. Consequently, there is no restriction on /i/i^ and /i can be 
taken arbitrary large. In the limit when /z — >■ oo one obtains exactly the 
same algorithm introduced in [21] (see also [IH]). In particular, one has 
Pmh,v = Pruh'^J^i and all that one needs to do is to solve for q = {I — Pm^)v, 
for which an explicit evolution equation is presented in [21] . 

Next, our aim is to give further conditions on h and n which guarantee 
that the difference between the approximating solution v and the reference 
solution u converges to zero as t — )■ oo. To do this we consider the time 
evolution of w = u — v. Since 

B{u, u) — B{v, v) = B{u, w) + B{w, v) = B{u, w) + B{w, u) — B{w, w) 

and 

Ih{u) ~ Ih{v) = Ih{w) 

then subtracting equation (125]) from equation (IT^ yields 

— + vAw + B(u, w) + B{w, u) - B{w, w) = -^iPJh{w). (32) 

This equation serves as the starting point for the proofs of Theorem [T] and 
Theorem [2] in the proceeding two sections. 

3 No-slip Dirichlet Boundary Conditions Case 

In this section we prove Theorem [TJ We first recall the following generalized 
Gronwall inequality proved in Jones and Titi [21] , see also [17] . 

Lemma 1 (Uniform Gronwall Inequality). Let T > be fixed. Suppose 

dY /■*+'^ 

— — h a{t)Y < 0, where limsup / a{s)ds > 7 > 0. 

dt t— >oo Jt 
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Then Y{t) — )■ exponentially, as t —)■ oo. 

We now state and prove a lemma leading to our main result. 

Proposition 1. Let Q be an open, bounded and connected set in H? with 
C^ boundary, and let u be a solution of the incompressible two-dimensional 
Navier-Stokes equations / fTPj) on Q with no-slip Dirichlet boundary condi- 
tions. Let V be the approximating solution given by equations (d^. Then 
|ti — f I — !■ 0, as t -^ oo, provided hcqH^ < v and fi > bc^G'^vXi. 

Proof. Let w = u — v. Then w satisfies equation ( 1521) stated above. Taking 
the inner product with w we obtain 

^-irkr + '^Ik'f + iB{w,u),w) = -fi{Ih{w),w) 
= fi{w - Ih{w),w) - fi\w\'^ 
< i^\P^(w - h{w))f + ^\w\' - fi\w\' 

/iCo/l^ll 1,2 Ai| |2 . ^u ||2 /^i |2 

< \\w\\ \w\ < —\\w\\ \w\ . 

- 2 " " 2' ' - 2" " 2' ' 

Since ( 1T6|) implies 

\{B{w,u),w)\ < c||m|||w|||w|| < — H-ulPlw^ H — ll^lp, 

we obtain 

d c^ 

— \w?+ (u Ikflkp <0. 

dt ^ V ^ 

Denote 

a{t) = /i ||w|p. 

V 

Taking T = (z/Ai)~^ in Theorem |4] we have for t >to that 

rt+T 



\\vf < 2(1 + Ti^\i)i^G^ = 4:i^G^ 



Thus 

rt+T 



limsup / a{s)ds > -^ - 4:C^G^ > c'^G^ > 0, 

t->-oo Jt VAi 

and by Lemma [1] it follows that \w\ — )■ 0, exponentially, as t — )■ oo. D 
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Proof of Theorem [II The hypothesis of Proposition [T] require that 
ficoh'^ < V and /i > Sc^G^z/Ai. 

Therefore 

-J^ > ^ > ciG^Ai (33) 

h^ V 

where c\ = Scqc^. D 

4 Periodic Boundary Conditions Case 

In this section we prove Theorem [2l We begin with an elementary inequahty 
which will be be referred to in the sequel. 

Lemma 2. Let (f){r) = r — /3(1 + logr) where /3 > 0. Then 

min{0(r) : r > 1} > — /31og/3. 
Proof. Note first that 

6(1) = 1-/3 and lim (h(r) = oo. 

r^-oo 

The derivative 0'(r) = 1 — f3/r is zero if and only if r = /3. Therefore 

1-/3 if < /3 < 1 



min{0(r) : r > 1} — 1 „, o -r i o 
'-^^ ^ - -' I - (3 log (3 if 1 < /3. 

Observe that over the interval < /3 < 1 we have 1 — /3 > — /31og/3, which 
concludes our proof. D 

We now state and prove a lemma leading to the proof of Theorem [21 

Proposition 2. Let Q = [0,L]^, for some fixed L > 0. Let u be a solu- 
tion of the incompressible two-dimensional Navier-Stokes equations [T^) on 
Q equipped with periodic boundary conditions. Let v be the approximating 
solution given by equations ^2B\j, where Ih satisfies ^. Then \\u — v\\ ^ , 
as t -^ oo, provided ficoh'^ < u and fi > 3z/Ai(2clog2c'^/^ + 8clog(l + G))G. 
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Proof. The proof makes use of the orthogonahty properties fll4p and (1151) 
along with the Brezis-Gallouet inequahty |S] which may be written as 

||M||L^(n) <c||m|| 1 + log ' , (34) 

which will allow us to obtain sharper estimates than for the case of no-slip 
boundary conditions. 

Take the inner product of equation ( l32l) with Aw and and use the orthog- 
onality relations (fT4|) and (TT5|) to obtain 

1 d\\wP 

2dt ^ ^^^^^^ ^ ^^^^' ^^' ^^^ ~ ^^(^h{w), Aw). 

Using ( 1M|) and the hypothesis fiCoh"^ < z/ we have 

\{B{w,w),Au)\ < c\\wfh + \og^^l^\\Aul 

and 

-~n{Ih{w), Aw) = fi{w - Ih{w), Aw) - fi\\w\\^ 

< fl\Pa{w - Ih{w))\\Aw\ -/i||w|P 

/i Cq/i 11 112 , '^i . |2 II ||2 / ^1 4 |2 A'li ||2 

< \\w\\ H — \Aw\ — ii\\w\\ < —\Aw\ \\w\\ . 



Therefore, 



at V I Ai\\w 



12 ( /^lll^ 

u,!. \ V, /vi||uy|| 

or 



|2 



(i||ii;p / \Aw\^ o I 4 ifi I^^P 



Now setting 



V Ai w r I AikuMJ / 



n 2c\Au\ , UwP 

p = ^ — and r = -;— ^^ — -^ 



in Lemma [21 and noting that r > 1, by Poincare's inequality ( ITTl) . we obtain 

d\\w\\^ r ^ , . ,, 2c|Au| 



+ /i-2c|AM|log ' \\\wf <0. 



dt 
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By ([22D we estimate 

, 2c\Au\ 

where 

J = C3 + C4log(l + G'), (35) 

C3 = 2c log 2c'^/^ and C4 = 8c. It follows that 



d\\w 



|2 



dt 
Furthermore, Young's inequality 



+ {/i- J|Au|}||^f <0. 



implies 



Denote 



J\M<f^\Au^'^l 



d\\wf If ^^ . ,2ln ,,2 



dt 2V fi 



J2 
Taking T = (z/Ai)~^ in Theorem H] we have for t > to that 



rt+T 



rt+i 

/ |Am|2<2(1+Tz/Ai)z/AiG2 = 4z/AiG2. 

Thus, 

limsup r^a{s)ds > -^ - ^J^G^ = -JG > 0, 
t^oo Jt 2z/Ai yU. 6 

and consequently \\w\\ — )■ exponentially, as t — )■ 00. D 

Proposition 3. Let Q = [0,L]^, /or some fixed L > 0. Let u be a solution 
of the incompressible two-dimensional Navier-Stokes equations ( fTPI) on Q, 
equipped with periodic boundary conditions. Let v be the approximating so- 
lution given by equations ^23^) where Ih satisfies ^. Then ||u — f || — )■ 0, as 
t ^ 00, provided ficoh'^ < v and [i > 3z/Ai(2clog2c^/^ + 8clog(l + G))G. 
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Proof. The proof is the same as the proof of Proposition [2] except that the 
estimate for —fi{Ih{w),Aw) has to be modified as 

-H{lh{w), Aw) = fi{w - Ih{w), Aw) - /i||tf IP 

< nlw — //i(w)||Au7| — /illwll^ 

4z/ 4 2 

Then, since —fi < — /i/2 the rest of the proof follows without change. D 

Proof of Theorem [H The hypothesis of Proposition |2] or Proposition [3] re- 
quire that 

fiCoh^ < V and /x > ?>v\\JG. 

Therefore, 

i_ > ^ > c2AiG(l + log(l + G)), (36) 

where C2 = 3max{c3, C4}. D 



5 Conclusions 

As shown in this paper, the algorithm given by ([2]), for constructing vit) 
from the observations Ih{u(t)), yields an approximation for u{t) such that 

\\u{t) — f(t) 11^2(5^) — )■ exponentially, ast — > 00, (37) 

provided the observations have fine enough spatial resolution. This result 
has the following consequence. To accurately predict u{t) for time T into the 
future it is sufficient to have observational data Ih{u{t)) accumulated over 
an interval of time linearly proportional to T in the immediate past. 

In particular, suppose it is desired to predict u{t) with accuracy e > on 
the interval [ti, ti + T*], where ti is the present time and T* > determines 
how far into the future to make the prediction. Let h be small enough and 
fi large enough so that Theorem [1] implies (1371) . Thus, there is a > and a 
constant C > such that 

\\u{t) - v{t)\\L2^n) < Ce""* for all t > 0. 

Now use v{ti) as the initial condition from which to make a future prediction. 
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Let w be a solution to (IT^ with initial condition w(ti) = v{ti). Known 
results on continuous dependence on initial conditions, see, for example, [H], 
[19] . [26] or [27], imply there is /3 > such that 

\\w{t) - u{t)\\L2^n) < \\w{h) - M(ti)||L2(n)e^(*-*i) for t > h. 

Therefore 

\\w{t) - u{t)\\ < Ce-''''+^^ < e for te[ti,ti+T] 

provided ati > f3T + ln(C/e). Thus w(t) predicts u(t) with accuracy e on 
the interval [ti, ti + T]. 

Work is currently underway to numerically test Theorem [2] in the case of 
determining finite volume elements and nodes. Of particular focus is how to 
tune the parameter /i. If fi is very large the effects of "spill over" into the 
fine scales may become significant, whereas if /i is small convergence of the 
approximate solution may be slow or not happen at all. Numerical simula- 
tions performed by Gesho [T8] confirm that the continuous data assimilation 
algorithm given by equation (^ directly works, without additional filtering, 
for observational measurements at a discrete set of nodal points, where Ih is 
given by dH]). As with previous computational work (cf. [I9],[2l] and [25] ) 
the approximating solution v{t) converges to the reference solution u{t) under 
much less stringent conditions than required by our theory. 

The main advantage of introducing a control term that forces the ap- 
proximate solution toward the reference solution is that we can rely on the 
viscous dissipation, already present in the dynamics, to filter the observa- 
tional data (that is, to suppress the spatial oscillations, i.e. the "spill over" 
into the fine scales, that are generated by the coarse-mesh stabilizing term 
filh{v)). In addition to working for a general class of interpolant observables 
this technique also allows processing of observational data which contains 
stochastic noise. In particular, the same algorithm can be used to obtain an 
approximation v{t) that converges (in some sense) to the reference solution 
u{t), to within an error of the order of /i times the variance of the noise in 
the measurements. This work j3] is in progress. 
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